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Over the grant period from October 1, 1992 to June 30, 1993 we performed the 

following investigations dealing with the generation of waves by contaminant ions. 

(1) Contaminant Ions Distribution Function: 

The properties of velocity distribution function of the contaminant ions were 
studied. Such ions are produced by charge exchange between ionospheric 0 + 
ions and the contaminant molecules and atoms, such as water (H 2 O). We 
followed the trajectories of a large number of H 20 + ions in a magnetic field. 
The initial velocity of the ions was chosen to be 8 km/s, which is the typical 
orbital velocity at the shuttle altitudes. The molecules and atoms released from 
the shuttle move with the orbital velocity with respect to the ionospheric 
plasma. We find that it takes about one ion-cyclotron period for the 
distribution function to transform from a beam (Figure la) to a ring distribution 
(Figure lc). 

(2) Linear Instability Analysis: 

We performed linear instability analysis to examine the possibility of wave 
excitation by an ion beam streaming perpendicular to the ambient magnetic 
field in an ionospheric type of plasma. The frequency spectrum over which 
growth of waves is likely to occur was determined parametrically depending on 
the relative density of the ion beam. We found growth can occur for 
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frequencies which are a sizable fraction of the ion-plasma frequency which can 

be as high as 50 kHz in the ionosphere. Figure 2a shows the wave frequency 
normalized to the ambient ion-plasma frequency {(Opo) as a function of 
no rmaliz ed wave number k x in units of cOpo / V Q where V Q is the orbital 
velocity. This figure is for Y = / Q e = 3, where (o pe and Q e are the 

electron plasma and cyclotron frequencies, respectively. The result shown is 
for propagation in a direction (0= 0) perpendicular to the magnetic field. 
Calculations have been performed for an appropriate range of 6. The 
corresponding growth rates (Im(a>)) of the waves are shown in Figure 2b. 
The growth rates are normalized to (Opo- The different curves in Figures 2a 

and 2b are for the relative densities of the water ions with respect to the 
ambient plasma density. It is seen that for a low beam density, the waves with 
the maximum growth rate have frequencies which are about 33% of the ion- 
plasma frequencies. For an ambient ion density of <10^ cm -3 , 
/ = 03 po / 2 tc < 50 kHz, making it possible for the excited waves to reach 

up to 17 kHz. It is also worth noting that the waves with maximum growth 
rates have nearly the maximum frequencies in curves for Re(<y) versus k x 
(Figure 2a); the maximum growth rate occurs near the peak of such curves. 

(3) 2 1/2-Dimensional PIC Code: 

We developed a 2 1/2 -dimensional PIC simulation code to study the nonlinear 
behavior of the waves excited by beam and ring types of ion distributions. The 
code treats both electrons and ions as particles. Since lower-hybrid waves may 
be involved, it is essential that both the particle species be treated kinetically. 
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(4) Nonlinear Studies of the Waves: 

We have performed initial simulation runs showing the excitation and 
saturation of the excited waves. In order to reach firm conclusions about the 
nonlinear behavior of the waves, a systematic set of runs need to be performed 
by varying the parameters of the problem such as the relative densities of the 
beam and the ambient plasma and the plasma temperatures. We are in the 
process of performing such simulations. 
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Figure 1: 
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Phase-space plots of ions produced by charge exchange, (a) At a very 
early time t = 1 iTj 1 , we see a beam-arc type of distribution. As more 
and more ions are produced, the beam-arc evolves into a ring 
distribution as shown in (b) and (c) for r = 5 and 12 CT e , 

respectively. 
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Figure 2: (a) Real part of the wave frequency (Re(fi>)) as a function of k x , the 

normalized wave number measured in units of (o^l V Q , where co^ 
is the ambient ion plasma frequency and V Q is the orbital velocity 
which is the initial velocity of water ions as they are released in the 
ionosphere. The curves shown are for different values of the 
normalized water ion densities, (b) The wave growth rates for the 
same parameter as in Figure 2a. The growth rate (Im(<y)) is given in 
units of (Opo- Note that the growth rate maximizes for a relative 

density °c = 0.5, i.e., when water and oxygen ion densities are the 
same. 





